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ABSTRACT 

We investigate numerically the contribution to the cosmic gamma-ray background 
from cosmic-rays ions and electrons accelerated at inter-galactic shocks associated with 
cosmological structure formation. We show that the kinetic energy of accretion flows in 
the low-red-shift inter-galactic medium is thermalized primarily through moderately 
strong shocks, which allow for an efficient conversion of shock ram pressure into cosmic- 
ray pressure. Cosmic-rays accelerated at these shocks produce a diffuse gamma-ray flux 
which is dominated by inverse Compton emission from electrons scattering off cosmic 
microwave background photons. Decay of neutral 7r-mesons generated in p-p inelastic 
collisions of the ionic cosmic-ray component with the thermal gas contribute about 30 
% of the computed emission. Based on experimental upper limits on the photon flux 
above 100 MeV from nearby clusters we constrain the efficiency of conversion of shock 
ram pressure into relativistic CR electrons to < 1%. Thus, we find that cosmic-rays 
of cosmological origin can generate an overall significant fraction of order 20 % and 
no more than 30 % of the measured gamma-ray background. 

Key words: acceleration of particles — cosmology: large-scale structure of universe 
— galaxies: clusters: general — gamma rays: theory — methods: numerical — radia- 
tion mechanism: non-thermal — shock waves 



1 INTRODUCTION 

After the re-ionization epoch, the inter-galactic medium is heated up to temperatures T ^ 10 4 K by energy input from quasars 
and star forming regions. As the wavelength of primordial density perturbations entering the non-linear stage, Xml(z), grows 
longer in size, the temperature of the low red-shift gas is further raised by shock waves. These convert the specific kinetic 
energy associated with accretion flows, of order \% L (z)H 2 (z) - with H(z) the Hubble parameter and z the cosmological 
red-shift - into thermal energy corresponding to 7-8 keV per baryon in massive clusters of galaxies and below ~ 1 keV in the 
diffuse inter-galactic medium (Cen & Ostriker 1999). The importance of this process in cosmic history is emphasized by the 
fact that on the one hand the large scale structure is most prominent in thermal X-ray radiation by hot intra-cluster gas. On 
the other hand numerical simulations predict that most of the baryons, which are observed at hi gh red-shift as Ly-q a bsorbers, 
end up shock-heated to temperatures in the range 10 s — 10 7 K in the contemporary universe (Cen & Ostriker 1999). 

Shock waves in cosmic environment are collision- less meaning that the dissipation process occurs via random electro- 
magnetic fields in the turbulent post-shock flow. In addition to a thermal population of hot gas, they generate a supra-thermal 
distribution of high energy particles, generally refer red to as cosmic- r ays (CR). Shock ac c eleration of CR electrons up to TeV 
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energ ies is directly obs erved in supernova remnants ( Allen et al. 1991] ; [Koyama et al. 1997 ; Tanimori et al. 1998 
2000, ; [Allen et al. 200 1| ) and the bulk of the Galactic CR ion population up to 10 — 10 ° cV is also thought to originate there 



The same acceleration mechanism is expected at work at large scale structure sh ocks. Indeed, we do observe non-thermal 
radiation in several massive clusters o f gala xies both at radio (e.g., Feretti 1999, and references therein) and hard X-ray 
frequencies ( Fusco-Femiano et al. 1999, 2000| , and references therein). At least in some models CR shock acceleration of both 



ions and electrons seems an essential ingredient to explain different features of the observed radio emission (e.g.. 



Miniati 



et al. 2001a). In addition, the possibility that CRs be accelerated at inter-galactic shocks has raised the iss ue that CR ions 



due to their long lifetime ag ainst energy losses and to their efficient confinement by magnetic irregularities ( Volk et al. 199£ ; 
Colafrancesco & Blasi 199£), may accumulate inside forming structures and store a significant fraction of the total pressure 
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there (Miniati 200C; Miniati et al. 2001b). Beside cosmic sh ocks, the energetics characterizing oth er CR sources such as 
termination shocks driven by galactic winds (Volk et al. 1996) and radio galaxies (EnBlin et al. 1997) are also noteworthy. 

Along these lines of observational and theoretical progress, the idea has recently been put forth that CR electrons 
accelerated at inter-galactic shocks may generat e the observed 7-ray b ackground (CGB) by way of IC (IC) scattering off 
cosmic microwave background (CMB) photons (Loeb & Waxman 200C). The CGB was measured by EGRET o nboard the 
Compt on G amma Ray Observatory. It consists of a differential flux, e J e ~ 5.4 (e /keV) -0 ' 1 keV cm -2 s _1 sr -1 (Sreekumai 
et al. 1 )98) throughout an energy range extending fr om about 100 MeV to 100 GeV. The integrated photon flux above 100 
MeV corresponds to 1.45 x!0~ 5 ph cm~ 2 s -1 sr -1 ( Brcckumar et al. 1995), in agreement with previously measures carried 
out by the SAS 2 satellite ^Thompson & Fichtel 198 



Its origin is most intriguing. The contrib ution from identified EGRET 
blazars to the integrated photon flux above 100 MeV amounts to lxl0~ 6 ph cm -2 s _1 sr -1 ( Mukherjee fc Chiang 1999 ); 
that is 



s sr 

7 % of the measured CGB flux in corresponding units. A number of authors have suggested that the contribution 
from an unresolved, faint end of the blazar 7-ray distribution function could fill the gap between the observed flux and that 
which is already accounted for. And in fact, though with large uncertainties, some authors were able to validate the above 
conjecture by using a 7-ray luminosity fun ction derived from the blazar radio luminosity function and the assumption of a 



constant ratio of radio to 7-ray fluxes (e.g., Padovani et al. 1993; Stecker & Salamon 1996). However, the correlation between 



t he blazar 7-ray an d radio flux is only weakly established and the existing evidence might indeed be dominated by bias effects 
(Miicke et al. 1997). In fact, based solely on the observationally deri ved luminosity function of 7-ray loud AGNs that have 
been detected by the EGRET experiment, |Chiang fc Mukherjee (1998 ) found that only 25 % of the CGB can be accounted for 
by undetected 7-ray loud blazars. Alternatively, the idea was put forward that the CGB is generated by the decay of neutral 
7r-mesons, it — » 77, generate d by an all pervasive population of CR ions interacting via p-p collisions with the nuclei of the 
diffuse inter-galactic medium (Dar & Shaviv 1995) . However, it was later argued that the CR ions from the available sources 
of CRs in the IGM, namel y shocks, Active Galactic Nuclei and normal galaxie s, could provide at most a fraction of order of 
a % of the observed CGB (Berezinsky et al. 1997; Colafrancesco & Blasi 1998). 

In this paper we carry out a numerical study to investigate the contribution to the CGB from CRs accelerated at 
cosmological, inter-galactic shocks associated with structure formation. We show that the heating of the IGM is due primarily 
to moderately strong shocks, which allow for some fraction of the shock ram pressure to be converted into CR pressure. 
According to our results, the contribution to the CGB from cosmological CRs could indeed be significant . It is dominated 
by I C e mission from primary electrons which, however, we estimate to be smaller than previously found (Loeb & Waxman 
200C()7} n addition we find a 7 -ray flux from 7r°-decay that is not negligible and larger than, although of the same order as, 
the upper limits estimated by |Colafrancesco fc Blasi (1998 ). 

The paper is organized as follows: we describe the numerical model in 
summarize with Sfi 



^, present the results in ^ discuss them in and 



2 THE SIMULATION 

Our investigation is based on numerical simulations of large scale structure formation which include shock acceleration, 
transport and energy losses/gains of CR ions and electrons. In the following sections we describe the salient features of the 
employed numerical techniques. 



2.1 Cosmological Model 



The formation and evolution of th e large scale stru cture is computed by means of an Eulerian, grid based Total- Variation- 
Diminishing hydro+N-body code ( Ryu et al. 1993[ ). The code is capable of sharply capturing both very strong and weak 



shocks while being computationally relatively inexpensive. This feature is of primary importance in this kind of study since 

cosmic shocks have direct impact on the CR properties. 

We focu s on a single, cur rently favored ACDM "concordance model" (Ostriker & Steinhardt 1995). The assumed model 
is flat (e.g., Jaffe et al. 2001) with a total mass density £l m = 0.3 and a va cuum energy de nsity 57a = 1 — fim = 0.7. The 
normalized Hubble constant is taken to be h = Ho/100 km s _1 Mpc -1 = 0. 67 (Frecdman 2000 ) and the baryonic mass density, 
fli, = 0.04. This value, after Big Bang nucleosynthesis calculations (e.g., Olive et al. 200C)is only marginally c onsistent (at 
the 2 a level) with the best fit to cosmic microwave background (CMB) measurements (e.g., laffe et al. 2001 ). The initial 
density perturbations are generated as a Gaussian random field with a power spectrum characterized by a spectral index 
n s = 1 and "cluster-normalization" as = 0.9. The initial velocity field is computed through the Zel'dovich approximation. 
The computational box is set to a comoving size L = 50 /i _1 Mpc. The dark matter component is described by 256 3 particles 
whereas the gas component is evolved on a comoving grid of 512 3 zones. Thus each numerical cell measures about 100 h~ 1 kpc 
(comoving) and each dark matter particle corresponds to 2 x 10 9 /i -1 M0. 



2.2 Cosmic-rays 



The CR dynamics is computed numerically through the code COSMOCR (Miniati 2001). The code follows simultaneously 
three CR populations, namely primary ions and electrons and secondary electrons/positrons (e ). Primary ions an d electrons 
are injected and accelerated at cosmic shocks according to the diffusive shock acceleration (DSA) mechanism (e.g.. 
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& Eichler 1987). Secondary e are produced in p-p inelastic collisions of CR ions with thermal IGM nuclei. In addition the 
code accounts for spatial tr anspo rt an d ener gy losses/gains undergone by each CR component. Before describing the numerical 
technique in some detail (§ |2.2.l| and § fc.2.2| ) we address some of the physical assumptions underlying our model. 

In general CRs are treated as passive quantities, meaning that at this stage their dynamical role is completely neglected 
both on the shock structure (test particle limit) and the gas dynamics. Recently developed kinetic models of DSA theory 
certai nly allow for solutions i n which shock acceleration is highly efficient and the CRs strongly affect the shock structure 
(e.g., Malkov & Drury 2001). In these solutions, however, most of the shock associated ram pressure is converted into CR 
pressure which becomes the dominant dynamic component. Alth ough it cannot be completely ruled out that such a situation 
occurs in the cosmic environmen t, both dynamical consideration s (Markevitch & Vikhlinin 1997; Horner et al. 1999| ; Nevalaincn 
et al. 2)0C; Roussel et al. 2000| Miralda-Escude Babul 1995; Wu"2000[ ) as well as existing upper limits on 7-ray emission 
from nearby clusters ( |Blasi 1999] ; [Miniati ct al. 2001b ), suggest that the ICM pressure is dominated by a thermal component, 
with room for a still significant CR pressure component at the level of a few tens of percent. The latter is likely to affect the 
gas-dynamics of the intra-cluster (and intergalactic medium) and will be considered in future work. 

Also, we implicitly assume the existence of a diffuse background magnetic field which is essential not only for the 
acceleration of the CRs but also for their subsequent confinement within cosmic structures. In support of this assu mption 
we noti ce t hat inter-galactic magnetic fields are indeed observed in the media of clusters of galaxies at the fiG level (Clarke 
et al. 2)01). Observations are much more difficult for the case of filame nts and therefore the evidence there is more sparse. 
However, the cas es of rotation meas ure from the Coma super-clusters (Kim et al. 198£) as well as the recent de tection of 
a ra dio filament (Bagchi et al. 2002) suggest that indeed sub-/iG magnetic fields might be common in the IGM (Kronberg 
1994|)~1 



2.2.1 Cosmic-ray Ions 



The first step in simulating CRs consist in detecting shocks where their injection and acceleration takes place. Shocks are 
identified when a converging flow (V ■ v < 0) experiences a pressure jump AP/P above a given threshold. For the present 
simulation, such threshold was set to the pressure jump across a shock with Mach number M=1.5. Successful detection of 
shocks requires a scheme that reproduces them as sharp transition in the numerical solution. This is especially important 
because relatively weak shocks, which tend to be harder to capture, cannot be ne glected. In fact, the y turn out to process a 



substantial fraction (but not all) of the kinetic energy associated to cosmic flows (Miniati et al. 2000, §3.1) 



The simulated CR ions are supposed to be injected in the DSA mechanism fro m the thermal IGM. Here the injection 
process is modeled after the thermal leakage prescription (e.g., Kang & Jones 1995). We assume that, upon shock passage, 



the gas thermalizes to a Maxwellian distribution defined by a post-shock temperature T s h oc k. Thermal ions in the high energy 
end of such a distribution mig ht be fast enough to "leak" back upstream of the shock against t rapping by the plasma waves 
that moderate the shock itself (Ellison & Eichler 1984; Kang & Jones 1995; Malkov & Volk 1995). These ions are thus injected 
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from the the rmal po ol in the DSA mecha nism and quickly gain energy much above the thermal average (e.g 
Eichler 1987). As in Miniati et al. (20011), based on cluster/group properties, we it found appropriate to follow CR ions up 
t o momenta p ma x — 10 D GeV/c. In fact, CR ions up to energies of 10 



10 1 



GeV/c can be accelerated at cosmic shocks 
([Kang ct al. 1996), althoug h their confinement within cosmic structures becomes difficult beyond 10 r GeV (Volk et al. 1996; 
Colafrancesco & Blasi 1998 ). The DSA process takes place on time-scales much shorter than the simulation time-step and, from 
a numerical point of view, it is treated as "instantaneous" . Thus, the injected particles are redistributed to a supra-thermal 
populatio n in accord w ith the shock properties. In the test-particle limit adopted here, this is a power law in momentum given 
by (e.g., iDrury 19831) 



/(p)shock = f(Pinj) Maxwell 



P 

Pin j 



(i) 



extending from p%nj to p m ax- The log-slope of the distribution is q = dlnf/dlnp = 3(7 9as + 1)/[2(1 - M~ 2 )] = 4/(1 - M~ 2 ) 
for 7 9 a S = 5/3 where M is the shock Mach number and 7 9 as is the gas adiabatic index. The Maxwellian and the power- law 
components that make up the full distribution function join smoothly roughly at the injection momentum, Pi n f, that is the 
momentum threshold above which ions are energetic enough to escape upstream the shock as mentioned above. The choice 
of pinj sets the normalization in the expression (111) . Here, as in previou s calculations, the momentum threshold is assumed as 
large as a few times the peak thermal value (e.g., Kang & Jones 199E), namely 



Pinj 



= Cl 2y/ mpfcsTshoc 



(2) 



where m v is the proton mass, /cb is the Boltzmann constant. In this study we adopt a parameter value ci ^ 2.5, on the 
relatively low-efficiency side of canonical values. The sign, >, indicates that when for very strong shocks (M ^ 10) the 
produced CR pressure is larger than allowed by the test particle approximation, we additionally impose an upper limit on it 
by demanding that it is ^ 40% of the shock ram pressure. This choice is admittedly artificial, although i t does resemble the 
effect of a CR back-reaction on the shock structure which suppresses the injection at the sub-shock (e.g 



Bcrczhko & Ellison 



199E). More importantly, however, it allows us to control how the adopted parameter affects the final result. 

In order to follow the further evolution of the injected CRs we divide momentum space into N p logarithmically equidistant 
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intervals (momentum bins) and for ea c h point mesh of the spatial grid, xj, we define the following piecewise power-law 
distribution function (Jones et al. 1999; Miniati 2001) 



/(x i) p)=/ i (x i )p-^^) J 



1 < Pj-i < P ^ Pj, 



(3) 



where pj... are the momentum bins' extrema. Spatial propagation and energy losses/gains of the accelerated CRs are then 
followed by solving numer ically a "kinet ic" equation written in comoving coordinates and integrated over the aforementioned 
momentum bins, namely ( Miniati 200 1|) 



dn(xj,p t ) 
dt 



-V- 



[un(xj,pi)]+ [b(p) 4ty P 2 /(x j ,p)]^_ i +Q(x s ,pi). 



(4) 



( j^) describes 

\dT J tot 



Here n(x,pi) = JJ* 47rp 2 /(xj,p) dp is the number density of CR in the i t h momentum bin, b(p) = 

adiabatic energy losses/gains as well as Coulomb collisions and p-p inelastic scattering and Qi is a source term, i(xj,pi), 
integrated over that momentum bin. At each time-step, gj (Xj) is determined self-consistently from n(xj,pj) and the required 
continuity of /(x 4 ,pj) at bins interfaces (Jones et al. 1999). Here we employed 4 momentum bins, although our tests indicate 
that the results are unchanged when N p is doubled. 



2.2.2 Semi-implicit COSMOCR 



Injection of thermal electrons in the acceleration mechanism is here simply modeled by assuming that the ratio between CR 



electrons and ions at relat ivistic energies is R e /V 
ran ge 1 x 10~ 2 - 5 x 10~ 2 flMiiller fc Tang 19871 



10 . Observationally, for Galactic CRs this ratio has been measured in the 



Vhiller & et al. 1995; Allen et al. 2001). The introduction of this parameter (see 



also Ellison et al. 200C) simplifies the treatment of this process, which is otherwise quite complex. In fact, the shock dynamics 
is presumably dominated by the proton component with which most of the mass, momentum and energy are associated. 
However, for an effective electron- wave interaction the resonant condition requires u> — kuvu = ±u>d, where ui and k are the 
wave frequency and wave- vector respectively, v is the particle speed, the symbol || indicates the component parallel to the 
ordered magnetic field, So, and uiL e ,i = eBo/m e ^c is the ionic/electronic (i, e) Larmor frequency. Thus, wave-modes generated 
by leaking protons, with k ~ OJu/v, are much lower than those with which a thermal electron can interact. Therefore, either 



elect r ons are pre-accelerated (e.g., McClemnts et al. 1997) or a different injection model is required (e.g., Levinson 199/ 



1994 [Bykov fc Uvarov 1999[ ). 

The dynamics of CR electrons is further complicated by the severe energy losses they suffer. For this reason the numerical 
treatment is far more delicate than for the CR ions as the distribution function steepens and cuts off above a certain momentum. 
In order to determine the evolution of this component more accurately, in addition to the electron bin number density described 
by eq. (|l|), we also follow the associated kinetic energy e, = 4-7T p 2 f(p) T(p) dp. Here T(p) = (7 — 1) m e c 2 is the particle 

kinetic energy and 7 = 1/ -y/l — (v/c) 2 is the Lorentz factor. The equation describing the evolution of e(xj,pi) is obtained in 
analogy to eq . (^), by integrating over the ith momentum bin a kinetic equation that has been multiplied by T(p). This reads 
([Miniati 200l|) 



de(xj,p t 
dt 



- V- [u e ( Xj , P i)] + [b(p) 4tt p 2 /(xj , p) T(p)] ™ 



Kp) 



47rep 3 /(xj,p) 
yj m e c 2 + p 2 



dp + S(xj,pi), 



(•») 



where S(xj,pi) = Air J p * p 2 T(p)i(xj,pi) and the third term on the right hand side, so written for numerical reasons, includes 
a combination of contributions from CR pressure work and sink terms. Thus the evolution of the CR electrons is described 
by eq. (Q) and (^) and now the slope of the distribution function at each grid point and momentum bin is determined self 
consistently by the values of n(xj,pi) and e(xj,p;). In addition, at each grid zone we also follow explicitly the upper momentum 

cut-off in the electron distribution function. 

For an accurate time integration of the system of eq. (^) and (JE)) the explicit formulation described in Miniati (2001) 
required a time-step, At, much shorter than the electron cooling time, r coo i. And in fact, in order to properly compute in a 
cosmological simulation (wher e At ~ 10 7 — 10 8 yr) t he electrons responsible for synchrotron and inverse Compton emission 
at the frequencies of interest, Miniati et al. (2001a) had to appropriately subcycle the integration of eq. (^) and (^) for 
several iterations. This procedure is not computationally convenient particularly at high red-shift. Furthermore, it becomes 
impractical when, in order to compute cosmic background radiation as in the present instance, the CR electron population is 
needed throughout most of the simulation. For this reason, for the electron population we have reformulated the momentum 
space portion (i.e., without the advection term) of the numerical scheme in a semi-implicit form. Such reformulation is quite 
straightforward an exercise and it mainly consists of expressing all the contribution to the fluxes in momentum space in terms 
of n(xj,pi) and e(xj,pi) at times t and t + At. The semi-implicit formulation is nominally less precise than the explicit one, 
although this will only affect the region near the momentum cut-off. Here, however, even though we explicitly follow the 
upper momentum cut-off, the major limitations of the numerical model are due to the sub-grid power-law description of the 
distribution function rather than to the semi-implicit formulation. Therefore, there is no significant loss of accuracy with our 
new approach. 

Given the finite simulation time-step, the time evolution of the very high and very low energy ends of the distribution 
function, characterized by a cooling time T coo i(p) < At, cannot be followed even with a semi-implicit scheme. We notice that 
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in the transport equation the diffusion term can be neglected away from shocks (Miniati 2001; Jones et al. 1999) and CRs are 



basically transported by the fluid (to which they are coupled through the magnetic field) at the flow speed v. Thus, because 
of the Courant condition, vAt ^ Ax, where Aa; is the mesh size, electrons such that r coo ((p) < At, will not propagate outside 
the grid zone in which they have been created. Therefore, in these energy ranges it is appropriate to take the steady state 
solution to eq. M) as 



P 2 /c(xj,p») dp = -4tt / — - 



dp 



p 2 i(x h p)dp. 



(6) 



We caution that for the case of shock acceleration, the above steady-state solution expresses the balance between a source 
and a loss term within the volume of a grid zone, but not at the acceleration region. In fact, the rate of energy gain at a 
shock must be higher than the loss rate up to the maximum energy of the accelerated particles. Physically, the expression in 
eq. (^) represents a cumulative population f c , i.e., a summation of all the individual populations of CR electrons within the 
grid zone xj that, after they emerge from the acceleration region, and as they are being advected away from the shock, start 
aging. Aging here refers to the modifications produced by energy losses. It is easy to infer from eq. (^) that if the source term 
is a power law with slope q s , the distribution f c will have an energy dependent slope 



9c (p) = 



91n/ c 



91n&(» 



(7) 



dlnp 9 hip 

For dominant synchrotron and inverse Compton losses, as in the case of very high energy electrons, eq. (0) gives q c ~ q s + 1, 
whereas when Coulomb collisions dominate q c ~ q s — 1,. Thus, effectively because of cooling and of a finite computational 
time-step, the numerical solution at a shock will be a broken power-law. At high energy the break will occur at momentum 



h 1fl 4 

Pbreak ~ 1U 



l iov; 



Ui 



I < MB ) 



MeV, 



(8) 



where Ub and Ucmb refer to the magnetic and CMB energy density respectively. At low energy, the break corresponds to 



i 

Pbreak 



(iov) ( 



io- 



MeV. 



(9) 



For the primary CR electrons we use 5 logarithmically equidistant bins, spanning momentum space from p m i„ = 15 MeV up 
to p 2 = 10 2 GeV. As for the case of CR ions, our tests show that the final results are not affected in any significant way by 
increasing the number of momentum bins. In addition, the distribution function above P2, computed according to the "steady 
state" solution [eq. (g)] for zones with non-null source terms, is recorded in an extra momentum bin ranging from p2 up to 
2 x 10 TeV. The latter is an appropri ate value for the maxim um energy of accelerated CR electrons, provided a magnetic field 
of order 0.1/iG throughout the IGM ([Loeb fc Waxman 200o|). 



Although not needed in the present study, we also have the capability of computing the sub-relativistic component below 
Pbreak- This component might be of interest, fo r example in order to estimate the emission of high energy radiation through 
non-thermal bremsstrahlung (e.g., Miniati 2002). Notice that in our model the energy associated to both CR ions and electrons 
is always fraction of the shock ram pressure. The small fraction of energy that is carried by the sub-relativistic components 
of both these CR populations, with time is transferred back to the gas in thermal form through Coulomb collisions. However, 
this same energy would have been converted into thermal energy through dissipation at the shock anyway, had DSA not be 
present. Therefore ; it is obvious th at the sub-relativistic populations created in our model do not produce any "thermodynamic 
catastrophe" (cf. Petrosian 2001, and references therein). 



2.2.3 Secondary Electrons and Positrons 

At each time-step we also compute the distribution of high energy 7-ray emitting secondary electrons and positrons (e 1 * 1 ) 
produced in hadronic collisions of CR ions with the nuclei of the intra-cluster gas, within the same momentum range as for 
the shock accelerated electrons. Secondary e^ 1 are generated in the decay of charged muons according to 



e +v e {y e ) +v ll {y v ) 



which in turn are produced in the following reactions 



P + P 
P + P 



>7V ± +X, 

K ± +X, 



p + v{v) 
* p ± + v(p) 

. ± 
+ TV + 7T 



(10) 

(11) 
(12) 

(13) 



In addition to p+p inelastic collisions, the above cascades are also triggered by the following interactions between CR ions and 
thermal nuclei: p+He, a+H and a+He. We include them by assuming a helium number fraction of 7.3% for the background 
gas and a ratio (H/He) ~ 15 at fixed energy -per-nucleon for the CRs (Vleyer et al. 1997). Further details on the source term 
for secondary are given in Miniati (2001). 
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Figure 1. Top: Two dimensional diagram showing the differential q lo g Mdlog Ti as a f unc ^ on of pre- shock temperature, Ti, and shock 
Mach number, M . AE t ^ is defined in eg. ([7^). It represents integrated thermal energy produced at cosmic shocks during cosmic history. 
It is shown in units of ke V per particle. 



3 RESULTS 

The shock strength, which we classify according to the shock Mach number, M, affects the CR produced background radiation 
in two respects: the normalization, which turns out higher for a flatter distribution function because it allows for a larger 
number of high energy, 7-ray emitting CRs; and the spectral index which, as shown below, has a simple dependence on the 
distribution function log-slope q. In addition, we notice that for a given injection rate, high Mach number shocks are more 
efficient in converting the shock ram pressure into CR ion pressure. As a first step, therefore, we investigate the strength of the 
shocks that were chiefly responsible for the heating of the inter-galactic medium and from which the primary CRs o riginate. 



The Mach number is given by the ratio of the shock speed to the local sound speed (Landau & Lifshitz 1987). When 
computing the shock Mach number we assume that the pre-shock gas temperature is not less than 10 4 K, the value it would 
quickly reach had photo-heating due to UV background light as well as radiation from the hot post-shock region been i ncluded 
For this rea son, the highest shock Mach numbers reported below do not exceed a few 100, contrary to the case in 



et al. (200C ) where such lower limit on the gas temperature had not been assumed. 
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3.1 Inter-galactic Shocks 

The amount of thermal energy produced at shock fronts throughout cosmic history is given by the expression 

/t(z=0) p 
dt / dSshockk B AT(nv) shock (14) 

where t is cosmic time, dS s hock is an element of shock surface and (nv) s hock and AT the gas flux and temperature jump across 
it. Thus, the quantity dS s hockkB^-T(nv) shock represents the amount of kinetic energy dissipated at the shock into thermal 
energy per particle and per unit time. In order to determine which shocks were primarily responsible for heating the IGM, 
in Fig. ^ we plot the differential quantity g g| jf§i^ f as a function of pre-shock temperature, Ti, and shock Mach number, 
M. It shows that in terms of parameter space (M, T), shock heating is confined to a narrow strip whose upper and lower 
edges are loci of constant velocity of order a few 10 3 km s _1 and 10 2 km s _1 respectively. A region where 3 < M < 10 and 
10 6 < T < a few x 10 7 further stands out where shock heating seems most "influential" . More specifically, from a quantitative 
analysis we find that thermal energy due to shock heating is contributed in the amount of 30% by weak shocks (M < 4), 45% 
by moderately strong shocks (4 ^ M ^ 10) and 25% by strong shocks (M > 10). Thus we expect most of the CRs accelerated 
at inter-galactic shocks to be described by flat distribution functions. 



3.2 The 7-ray Background 

The 7-ray flux in units of 'keV cm -2 s _1 sr -1 ' at a given photon energy, e, is computed as 
s rM _ _c /"""" e ~ T " 7 j[e{l + z),z] 

eJ(£) - £ 4^J [n m (i + z )3 + n A ]i/2 (i + zY dz (lb > 

where j(s, z) is the computational-box-averaged spectral emissivity in units 'photons cm -3 s -1 ' computed at each simulation 
red-shift, z, and at the appropriately blue-shifted photon energy e(1 + z). In addition, r 77 is an attenuation factor due to 
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Figure 2. Top: The totals-ray background flux produced by cosmological CRs (solid line). It includes the contributions from IC emission 
from shock accelerated CR electrons (dotted line), decay of n° generated in inelastic p-p collisions ( dashed line ) and IC emission from 



secondary e ± (dashed- dotted line). The experimental data and their error-bars (solid dots) are taken from Sreekumar et al. (199t ) 



photo-pair creation, 77 — > e . It can actually be ignored given that the energy, e, of the measured 7-ray photons is below 100 
GeV and that the contribution from the high red-shift universe is not significant (see below). Finally, c indicates the speed 
of light and z max an upper limit of integration. The actual value of the latter parameter is not too important as most of 
the 7-ray emission is produced in the low red-shift universe, with j(s, z) decreasing by about an order of magnitude per unit 
red-shift increase. The strong red-shift evolution of j(e, z) is related to the available am ount of kinetic energy to be dissipated 
at shocks, which is itself a strongly decreasing function of red-shift (Miniati et al. 2000). We consider two emission processes: 
IC of CR electrons scattering off CMB photons (leptonic), and decay of neutral pions, it — > 77 (hadronic), produced in p-p 
inelastic collisions of CR ions off nuclei of inter-galactic gas. In addition we also consider IC emission from secondary 1 



produced in the aforementioned p- p collisions. Our treatment of IC emission and 7r° production are detailed in 
(2001a| and "~ 



Miniati et al 



Miniati et al. (2001b) respectively. For IC emission we simply use a Thompson cross section, as Klein-Nishina 



corrections become relevant only for electron energies about an order of magnitude above the maximum considered here. The 
results are shown in Fig. ^ where we plot e 2 J(e) as a function of e. The total flux (solid line) corresponds to a constant value 
at the level of 0.2 keV cm~ 2 s _1 sr _1 throughout the spectrum. It is clearly dominated by IC emission from shock accelerated 
CR electrons (dotted line). A fraction of order 30% is also produced by 7r°-decay (dashed line), whereas IC emission from 
secondary e* (dotted-dashed line) turns out about an order of magnitude smaller, although it is of the same order as 7r°-decay 
and becomes larger th an it around 100 MeV. For comparison the observational data, taken from the EGRET experiment 
( sreekumar et al. 1995 ) , are also plotted with the associated error-bars (solid dots) . 

Remarkably, all three components produce the same flat spectrum, similar in shape to the observed one, except below ~ 
GeV where the ty° production cross section drops due to phase-space constraints. This result is a reflection of the fact that the 
CRs at the origin of the 7-ray radiation were produced in strong shocks and have a flat distribution function with a log-slope 
q ~ 4. Indeed, IC emissivity from a power-law distribution of CR electrons, f(p) oc p~ q ° , is of the form j(e) oc e - ^ -3 )/ 2 . 
Therefore, for a steady state distribution function of CR electrons accelerated at strong shocks (Al > 4) and subject to IC 
losses q e ~ 5 [see eq. (Q)], e 2 J(e) oc ej(e) oc e°. As for the hadronic component, the pion emissivity (in units of pions s~ 
cm -3 ) from a 



[see eq. (Q)], e 2 J(e) oc ej(e) oc e u 
CR ion distribution f(p) oc p" qp , is j^ T o(e 7l 



-(4q p -13)/3 



(cf. [Mannheim fc Schlickeiser 1994[ ). Thus, for 



q v ~ 4 (strong shocks), as before e J(e) oc e . Analogously, the secondary e 



a stead y state distr ibution against IC losses of the form f e ± (e e ± ) oc e ± 
& Schli 



(4<j p 



injected in the above hadronic processes reach 
for q p ~ 4 [cf. eq. (@) 



OC £ ± 



and 



Mannheim 



Jceiser 1994 . As the steady state distribution of primary e , this component also produces an IC emission j(e) oc e 



and a 7-ray background flux ej(e) oc e°. 

Although with a spectral shape that accords with observations, the computed 7-ray flux is only 15 % of the observed 
CGB. A number of factors, which we address below, might affect these results. Firstly, resolution effects could be important, 
particularly for the hadronic component, because neutral ir- mesons are produced in a two body process. In fact the emission 
from smaller structures is likely to be underestimated, although we do not expect the hadronic component to be dominated 
by the contribution of the low end of the mass function. This is corroborated by the fact that our present estimate for the 
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Figure 3. Top: The IC ■y-ray photon luminosity above 100 MeV from individual clusters as a function of the cluster X-ray luminosity 
weighted core temperature. The flux is produced by CR electrons accelerated at accretion shocks associated to groups/clusters of galaxies. 



flux produced by this process is alrea dy larger than, although of th e same order of magnitude as, the upper limit given by the 
analytical calculations carried out by Colafrancesco & Blasi (199£). In particular, between 1 GeV and 100 GeV the computed 

7% when combined with the contribution 
Colafrancesco & Blasi (1998) estimated a range of values 0.5 — 2 % (but they did not 



7-ray flux from 7r°-decay corresponds t o about 5 % of the measured value and to 
In comparison 



of IC emission from e 

include IC emission from secondary e^). Thus the error on the 7-ray flux as a result of hadronic processes should amount to 
a factor two or so. As for shock accelerated electrons, our tests indicate that coarse grid effects are not important and that 
they amount at most to a few tens of percent. 

More important is the issue related to the assumed acceleration efficiency, 0, for converting the shock ram pressure into 
CR pressure. In this calculations we adopted Qmax ^ 40% with only a small fraction R e / P ~ 10~ 2 of the generated CR pressure 
to be borne by electrons. Shock acceleration can be even more efficient than allowed here. However, in those cases we expect 
non-linear effects to suppress the injection rate with consequent reduction of the 7-ray emitting CR population whilst most 
of the CR pressure is carried by the upper end of a concave-up distribution of CR ions. In principle the parameter R e / P could 
be higher, allowing for a substantially larger contribution from the leptonic component. This point is particularly important 
because it has been re cently claimed that IC 7-ray flux from CR accelerated at inter-galactic shocks can make about all of 
the unexplained CGB (Loeb & Waxman 200C). Therefore, in the following we will try to assess the efficiency of conversion of 
shock ram pressure into CR electrons based on EGRET observations of nearby clusters. 



3.3 IC 7-ray Emission from Clusters of Galaxies 

For this purpose, we have computed the integrated photon 7-ray luminosity above 100 MeV, L 7 (> 100 MeV), produced by 
primary CR electrons accelerated at accretion shocks associated to groups/clusters of galaxies. L 7 was computed within a 
volume of r adius about 5 ft _1 Mpc around the cluster center, although most of the flux is produced roughly within a radius of 
3 fe~ 1 Mp c (Miniati 2002). C ollapsed objects w ere identified throug h a slightly modified version of the "spherical over-density" 
method (Lacey & Cole 1994), fully detailed in Miniati et al. (200C). However, because of the large spatial extension of the IC 
emission, there is a real risk of mis-association with any selected object of radiation emission produced in nearby but unrelated 
accretion shocks. To circumvent this problem we further refined our sample by rejecting the smaller of two or more objects 
that happen to be closer to each other than a distance of 7.5 Mpc. In fig. ^ we plot the computed 7-ray photon luminosity, 
L 7 (> 100 MeV), against T x , that is the thermal X-ray luminosity weighted cluster core temperature. The correlation between 
the plotted quantities is tight and the numerical data are well fit with a single power-law function. After performing a simple 
least-x 2 analysis we obtain the best-fit parameters as 



L 7 (> 100 MeV) 



8.7 x 10 4 



R, 



Iv 



10- 



(— v 



ph s 



(16) 
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Table 1. 



Object 


T x 
keV 


d L 
Mpc 


F 7 (>100 MeV) a 
ph cm -2 s _1 


References 


Virgo 


1.8 


15.5 


< 2.20 X 10" 8 


1,2,3 


Perseus 


6.3 


82.5 


< 3.73 x 10~ 8 


4 


Coma 


8.3 


104.9 


< 3.26 x 10~ 8 


5,6 



a From O. Reimer, priv. comm. 

Reference s - (1) Graham et al. 1999 ; (2) Fouquc ct al. 2001 (3) Bohringer et al. 1994 (4) 
2001; (6) [Baum ctal. 199?j ; 



Schwarz et al. 1992 



(5) Arnaud & Evrard 



The temperature scaling for L 7 obtained is close to what expected based on clusters scaling relations. In fact, in a steady state 
configuration appropriate here due to the short CR electrons cooling time, the radiation flux must be produced at the same 
rate at which kinetic energy is being supplied to the external accretion shocks. That implies F y tx pvi hock B? ~ T 2 ' 5 . The plot 
in fig. ^ is noteworthy by itself. In fact it shows that the 7-ray flux associated to IC emission from shock accelerated electrons 
is comparable with the flux from 7r°-decay which is usually regarded as the main source of 7-rays from clusters. Thus, for a 
correct interpretation of the 7-ray emission from clusters that planned 7-ray f acilities aim a t detecting, this component needs 
to be properly addressed. This task is being undertaken in a separate paper (Miniati 2002). 

In order to constrain the electron acceleration efficiency in fig. ^ we compare the computed 7-ray photon luminosity above 
100 MeV to upper limits s et by the EGR ET experiment on the 7-ray flux from nearby clusters of galaxies. Of the 58 cases 
that appear in the list by Reimer (1999), we found that the tightest constraints were those set by the experimental upper 
limits on the following clusters: Virgo, Perseus and Coma. The latter two have temperatures larger than those of objects found 
in the computational box and, therefore, their upper limits must be compared to values extrapolated from the the scaling 
relation given in eq. (|l6|). The flux upper limits were converted into luminosity limits based on the relation L 7 = 47rci 2 F 7 , 
where d,L is the luminosity distance taken, for each cluster, from the literature. In table |l| we summarize the assumed values 
for the cluster temperature, dist ance and 7-ray flux upper limits. The observational data adopted here, although numerically 
very close to those published in Reimer (1999), were recently re-obtained by Reimer after careful and improved analysis of 
numerous staked observations for each source (Reimer, priv. comm.). The plotted experimental upper limits (particularly 
those relative to Coma and Virgo clusters) require that the fraction of shock energy that goes into relativistic CR electrons 
be ^0.8% x [log (p max /m e c)/ log(4 x 10 7 )] where p ma x is the maximum energy of the accelerated electrons. For the adopted 
ion injection efficiency this roughly translates into R e / P ^ 2 x 10~ 2 . In any case, this allows us to set an upper limit on the 
computed 7-ray flux of about 0.35 keV cm -2 s _1 sr~ . Thus, according to our computation, cosmological CRs could overall 
contribute an importan t fraction of order ~ 25 % of the extrajjalactic cosmic 7-ray background as measured by EGRET 
( Sreekumar et al. 1998| ). Due to the scatter in the plot of fig. which we regard as "real" and attribute to the different 
shock histories of the simulated objects - it is plausible to relax the upper limit on the acceleration efficiency of CR electrons 
by another factor < 2. This is allowed, however, only if the 7-ray flux produced by 7r°-decay is negligible. And although this 
is an extreme scenario, it cannot be excluded because 7-rays from galaxy clusters/groups have never been measured. Thus in 
this latter case the IC 7-ray flux could account for up to ~ 30 % of the measured CGB. Notice that the upper limits derived 
here rely mainly on the assumption that we compute correctly the relative fraction of 7-ray flux from accretion shocks onto 
clusters and on smaller structures, where most of the emission actually originates. 

To summarize, the whole unexplained 7-ray background flux could be accounted for by assuming an acceleration efficiency 
for the electronic component ~ 3%, or higher by a factor ~ 7.5 than assumed in fig. ^, but that would result in fluxes from 
nearby clusters such as Virgo, Coma and Perseus well above th e experimental upper limits. It is instructive to compare in some 
more detail our approach and that of Loeb & Waxman (2000) as the differences go beyond the assumed acceleration efficiency 



for the electronic component. In fact, in Loeb & Waxman model, the population of emitting CR electrons is computed by 
postulating that a fraction £ e of the thermal energy in the shocked IGM is converted into a flat power-law of relativistic 
electrons. The thermal energy density of the shocked IGM is estimated as 1.5f s h,n p kBT, where f s h is the fraction of shocked 
baryons, n p the average proton n umber density and T the mass-weighted average gas temperature. For their quantitative 
estimates Loeb & Waxman (2000) assume f a h(kBT /IkeV) ~ 1 and £ e ~ 0.05. Notice that, unlike £ e , the efficiency parameter 



defined in our case is based on the shock ram pressure instead of the post-shock thermal energy density. For moderately 
strong and strong shocks, the two definitions differs by 12.5 % (it gets more complicated for weak s hocks but this case is no t 
relevant because weak shocks are inefficient accelerators). In particular, assuming £ e ~ 0.05 as in Loeb & Waxman (200C) 
means that 5.6 % of the ram pressure is converted into CR electrons. This corresponds to a factor ~ 15 — 19 larger than the 
value of 0.3-0.4 % adopted to compute the dotted line in fig. ^| which lies, however, only a factor 10 below Loeb & Waxman 
's estimate. This difference by a factor ~ 1.5 — 2 gro ws if one accounts fo r the following additional matter s: ef fectively in the 
simulation f s h ~ 2/3 and ksT/keV ~ 1/2 (see also Keshet et al. 2002). In addition, as pointed out in [3.1, only ~ 2/3 of 
the IGM thermal energy comes from shocks strong enough (i.e., M 4) to produce 7-ray emitting CR electrons. Thus, using 
the appropriate fudge factors Loeb & Waxman expression would produce a result smaller than the simulation by almost an 
order of magnitude (factor ~ 7). In part this is due to the fact that, by following the full evolution of the shock accelerated 
CR electrons from 10 MeV to 2 TeV, we are also able to account for a re-acceleration process. That is the re-energization 
of aged (by cooling) CR electrons that were produced in earlier shock events. As already pointed out, this effect increases 
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Figure 4. Top: Histogram of the integrated photon flux above 100 MeV against temperature of the IGM where the emission originates. 
We show the case for IC emission from shock accelerated CR electrons (dotted line), for emission from n° -decay (solid line) and for IC 
emission from secondary (dashed line). 



our computed CGB by a factor between 2 and 3. Part of the remaining discrepancy could be due to 7-ray flux emitted 
at high red-shift (z > 0.5) that Loeb & Waxman neglected, although according to our simulation that should not be too 
significant. Most likely, however, it reflects the fact that the actual conversion of shock energy into CR electrons cannot be 
simply estimated through the IGM average temperature. In any case, given its simple formulation that estimate is quite close 
to our results. 

Just prior to submission of this pap er, the work by Keshet et al. (2002) was brought to our attention. The authors further 
investigate the Loeb fc Waxman (2000| ) model utilizing smoothed particle hydrodynamic (SPH) numerical simulations. They 
detect the occurrence of shocks in individual SPH particles by measuring changes in their entropy every some time interval. 
When shocks are identified a fraction £ e of the generated post-shock thermal energy is converted into CR electron energy. 
The CR electrons are distributed in energy according to the test particle limit of DSA theory, as assumed here. Their results, 
particularly the level of CGB produced by shock accelerated electrons, is in good agreement with our estimated flux plotted 
in fig. [| However, there are some differences that are worth pointing out. In fact, for the same estimated 7-ray flux our 
simulation requires an electron acceleration efficiency smaller than theirs by about an order of magnitude (of course this 
includes already the different definitions of efficiency) . We find that most of the thermal energy in the low red-shift universe 
was produced by what we refer to as moderately strong shocks, that is shocks with Mach number in the ran ge 4 -10. Those 
authors mention some difficulty in tracking these shocks. If these are neglected according to our results in £3.1 one would 



be considering only about 1/3 of the total flux through shocks (although their estimate indicates more of a factor 1/2). In 
addition our simulation also includes a re-acceleration process which, as mentioned before, boosts the computed flux by a 
factor 2-3. These two factors, namely that related to moderately strong shocks and re-acceleration, are probably the most 
important in determining the differences between our and Keshet et al. results. 



3.4 Emission Isotropy 

In order to assess the spatial distribution of the 7-ray emissivity, in fig. ^ we show a histogram of the total integrated 
photon flux above 100 MeV, F 7 (> 100 MeV)= J 100McV J{e)de, versus the temperature of the inter-galactic medium where 
it originates for the three computed components. The diagram illustrates that the contributions from both 7r°-decay (solid 
line) and secondary e* (dashed line) are strongly concentrated within regions with temperature T ^ 10 7 keV, corresponding 
to groups and clusters of galaxies. On the contrary the IC emission by shock accelerated CR electrons (dotted line) is roughly 
equally produced in the shocked IGM with temperatures ranging from a few xlO 5 K up to several xlO 7 K. Notice, however, 
that because of strong radiative losses, these high energy electrons are almost exclusively found in "narrow" post-shock regions. 
Thus for this component (dotted curve), the temperature axis in fig. ^ is indicative of the gas temperature behind shocks, 
rather than anywhere throughout the diffuse IGM. In any case, the radiation flux produced by this component is expected 
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to be highly isotropic (Waxman & Loeb 200C) and, as it originates from the shocked diffuse IGM, if detected will allow an 
independent lower limit on the cosmic baryon density (Loeb & Waxman 2000). 

For all components, the 7-ra y surface brightness will turn out most prominent toward galaxy clusters even for the IC 
emission. Waxman & Loeb (2000) predict for the IC component fluctuations of the 7-ray background intensity of order 40 % 
on the sub-degree scale. This fluctuations should actually be even larger when the hadronic component that they neglected is 
accounted for. I n fact, the latt er should be apparent only toward massive clusters and should dominate the emission within 
the cluster core (Miniati 2002). 



4 DISCUSSION 



We have constrained the energy acceleration efficiency of relativistic electrons to < 1 x [\og(p rnax / m e c) / log(4 x 10 7 )] % 
based on EGRET experimental upper limits on the photon flux above 100 MeV from nearby clusters. We specify that 
the found efficiency limit is valid only insofar as the accelerated CR electrons are power-law distributions extending above 
Pmax ~ \/ 75MeV / (ecmb ) ni e c ~ 100 GeV, where (ecmb) is the average CMB photon energy. And, in fact, strictly speaking 
the found constraint actually regards the total number of CR electrons in this energy range. This upper limit on the electron 
acceleration efficiency plays a crucial role in the determination of the allowed contribution from cosmological CRs to the CGB. 
One could doubt that the EGRET measurements would apply for extended sources such as galaxy clusters accretion shocks, 
which measure several degree on the sky. However, the EGRET field of view is much larger than this and the exposure times 
drop by only a factor of two at about 18° from the instrument axis. Thus, had there been any extended flux above the upper 
limit set by the EGRET experiment they would have been detected (O. Reimer priv. comm.). 

The upper limit just obtained on electron acceleration efficiency strictly speaking concerns the acceleration at cluster 
shocks, whereas a significant part of the 7-ray emission is produced in inter-galactic shocks. However, since clusters of galaxies 
are the deepest potential wells in cosmic environment, their associated shocks are indeed the strongest. Therefore, if anything, 
the emissivity is expected to be highest there, unless these shocks are strongly modified by a back-reaction of CR pressure. 
In addition, unlike for massive clusters, shocks from accretion flows on filaments and smaller structures in general could be 
weakened by energy input into the IGM from supernovae, stellar winds and radio-galaxies. These non-gravitational processes 
are thought to raise the temperature of the IGM up to 10 7 K , as it seems required by the existence of an "entropy floor" in 
the scaling relation of groups of clusters (Ponman et al. 1999). In that ca se, several inter-galact ic shocks would be suppressed 
and the 7-ray background flux generated there substantially reduced ( Ibtani fc Inoue 2002] ). From our fig. ^ we infer a 
correction factor with res pect to our estimate of order ~ 1/3. However, the actual role of such feed-back processes is still 
highly controversial (e.g., Kravtsov & Yepes 2000, and references therein) and alternati ve scenarios in which the presence of 
cooling reduces dramatically the amount of required pre-heating are being investigated (Voigt & Bryan 2001, and references 
therein). 

According to our result in fig. |j| the the Gamma-ray Large Array Space Telescope should be able to measure the 7-ray 
flux due to IC emission from CR electrons accelerated at cluster accretion shocks (see also Miniati 2002). In a previous 



work we explored with the same numerical technique adopted here the role of cluster accretion s hocks as sites where the CR 
electrons that power radio relics (and not radio halos) are being produced (Miniati et al. 2001a). Our numerical model was 



able to account for several observed emission properties of radio relics (see also Enfilin et al. 1998; Roettiger et al. 1999, for 
direct observations related to this model). Interestingly the upper limit on the electron acceleration efficiency found here, 
corresponding to a parameter R e / P ~ 2 x 10 -2 , is in line with the values adopted in that work in order to reproduce the 
measured ra dio emission for a number of sources. This number is also in agreement with recent observations of supernova 
remnants by Allen et al. (2001), who find R e / P ~ 1/160. In any case, upcoming observations both in soft and hard 7-ray 
should allow us to put better constraints on such parameters for clusters accretion shocks. 

Thus, with these constraints on the acceleration parameters the contribution to the CGB from cosmological CRs is 
estimated at the level of 20 % and most li kely below 30 % of the measured flux. We point out that according to recent 
detailed calculations of [Strong et al. (2000| ) the diffuse continuum 7-ray emission from the Galaxy could be significantly 
higher than previously estimated, thus lowering the actual level of the diffuse extra-galactic CGB. Likely, the next generation 
of high sensitivity 7-ray observatories will be resolve part of the CGB into a faint population of AGNs. In any case, the flux 
contributed by cosmological CRs will perhaps turn out the only truly diffuse extra-galactic 7-ray background component. 



5 CONCLUSIONS 

We have investigated numerically the contribution to the CGB from cosmological CRs. We carried out a simulation of structure 
formation in the canonical ACDM universe and followed directly the acceleration, spatial transport and energy losses/gains of 
three different CR components: shock accelerated ions and electrons as well as secondary e^ 1 . Our results can be summarized 
as follows: 

• Cosmological CRs produce a significant fraction of order of 20 % and no more than 30 % of the measure CGB. The 
computed flux is dominated by a leptonic component (~ 70 %) and accounts for a non negligible contribution from hadronic 
processes (~ 30 %). 
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• A higher flux, although in principle admissible, would be in disagreement with experimental upper limits set by EGRET 
on 7-ray emission from nearby clusters of galaxies. 

• Based on the above considerations, we set an upper limit on the efficiency of conversion of shock ram pressure into 
energy of relativistic CR electrons of order ~ 1 x [log(p maa; /m e c)/log(4 x 10 7 )] %, where p max is the maximum energy of the 
accelerated electrons. 

• IC 7-ray emission from clusters of galaxies could be of the same order of magnitude as 7-ray emission from 7r°-decay. 
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